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Abstract. We study a generalization of the Wolf- Villain (WV) interface growth 
model based on a probabilistic growth rule. In the WV model, particles are 
randomly deposited onto a substrate and subsequently move to a position nearby 
where the binding is strongest. We introduce a growth probability which is 
proportional to a power of the number 7ij of bindings of the site i: pi <x n". 
Through extensively simulations, in (1 + l)-dimcnsions, we find three behavior 
depending of the v value: t) if v is small, a crossover from the Mullins-Hcring 
to the Edwards- Wilkinson (EW) universality class; it) for intermediate values of 
v, a crossover from the EW to the Kardar-Parisi-Zhang (KPZ) universality class; 
iii) and, finally, for large v values, the system is always in the KPZ class. In 
(2 + l)-dimcnsions, we obtain three different behaviors: i) a crossover from the 
Villain-Lai-Das Sarma to the EW universality class, for small v values; it) the 
EW class is always present, for intermediate v values; iii) a deviation from the 
EW class is observed, for large v values. 
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1. Introduction 

As the interface growth is ubiquitous in nature, it have constituted an important 
subject of research in the past decades [TJ [2]. In general, the thickness of 
these interfaces, also named roughness w, is a signature of the nonequilibrium 
growth conditions under which the interface has formed and their evolution can be 
characterized by the Family- Vicsek scaling relation [3] 

v>(L,t)=L°f(±^ , (1) 

where L is the linear size of the system and t is the evolution time. The scaling 
function f(x) behaves as 

fix) ~ { ^' f ° r X <<C 1 (2) 
w \ const, for x > 1 w 

where /3 — a/z. Then, the roughness w grows with t@ until saturates in a value w sat 
denominated saturation roughness, which behaves with the system size L as a power 
law according as L a . The saturation times t sat also grows with the system size as L z . 
The exponents a, /3 and z are known as roughness, growth and dynamic exponents, 
respectively. 

This scaling analysis lead to a considerable advance in the understanding of the 
roughening of interface growth and allow to define some universality classes which 
are related to a stochastic differential equation for the height h(x, t) of the system 
at position x and time t. Lai and Das Sarma [3] propose a general form for these 
equations 

dh 

— = v V 2 h+v 1 V i h+. . .+A (V/ l ) 2 +A 1 (V 2 h) 2 +A 2 V-(V/ l ) 3 +. . .+7] , (3) 

where the parameters Vi are related to the linear terms, A^ to the nonlinear ones and 
X] is a white noise. 

In the hydrodynamic limit [L — > oo and t — > oo), the Ao term dominates and the 
asymptotic behavior is governed by the KPZ universality class [5]. This equation is 
the simplest nonlinear differential equation to describe a kinetic growth process and it 
has exact solution only in (1 + l)-dimensions where the scaling exponents are a = 1/3 
and z — 3/2. 

For models where the up-down symmetry is present, Ao = [I] and the vq term 
of equation Q dominates. In the asymptotic limit, we recuperate the linear equation 
proposed by Edwards and Wilkinson [3] to studied a sedimentation process. This 
equation has an exact solution and the scaling exponents are given by a = (2 — d)/2, 
where d is the dimension of the substrate, and z = 2. 

However, for short and intermediate times, the process can present a behavior 
which is a portrait of the dominant term before of the hydrodynamic limit. An 
example is the linear equation proposed by Wolf and Villain 7] and by Das Sarma 
and Tamborenea [H] to introduce surface diffusion where the term i>i dominates for 
short times. Equation ^ just with these term is linear and the exact solution gives 
a = (4 — d)/2 and z = 4. In table 1 we show the growth and the roughness exponents 
associated to the main terms of equation in d = 1 and d = 2. 

The determination of the asymptotic universality class of a growth model may be 
a difficult task due to the extremely long time of the initial behavior that can mask the 
crossover to another one [9^ ■ To investigate this problem various techniques have been 
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Table 1. The growth /3 and the roughness a exponents associated to the main 
terms of equation (2) (left column) for d = 1 (center) and d = 2 (right) QQ. 





(3 a 


f3 a 


V 2 h 
(Vhf 

V 4 h 
V 2 (Vh) 2 
V- (V/i) 3 


1/4 1/2 
1/3 1/2 
3/8 3/2 
1/3 1 
3/10 3/4 


(log) (log) 
w 0.24 w 0.4 
1/4 1 
1/5 2/3 
1/6 1/2 



developed, as examples we can cite the investigation of the universality class using 
the surface diffusion currents [101 lllj , the application of the noise reduction technique 
to models with limited mobility [12], the investigation of finite size effects in the DT 
and WV models |13j . the application of the scaling transformation to the stochastic 
equation of the WV model [HI [T3] . 

In this article, we propose a probabilistic interface growth model that is a 
generalization of the WV model where an adatom diffuses to maximize its coordination 
number. In our model, the incoming particle search for a site where the number of 
bonds is higher and it has a probability to be incorporate in this site. The probability 
depends of the number of bonds that the particle will have and it has a parameter v 
so that, for v = 0, the WV model is recuperated. For low values of z/, the crossover 
to the asymptotic universality class occurs early and, for high values, the non-linear 
KPZ universality class appears. A similar strategy was usaed in the investigation 
of the on-lattice Eden model [16], where the authors show that the introduction of 
such probability growth leads to the growth of on-lattice Eden clusters whithout the 
undesirable anisotropy effects. This strategy was initially proposed to generate an 
isotropic cluster of the on-lattice diffusion limited aggregation model [T7]- However, 
a later work [TH] shows that in the noiseless limit, instead of isotropic patterns, a 45° 
(30°) rotation in the anisotropy directions of the clusters grown on squre (triangular) 
lattices was observed. The model that we study is presented in the next section and 
the results in section 3. The last section shows some conclusions. 



2. Probabilistic Growth Rule 



We study a model where the particles are randomly deposited on a d-dimensional 
lattice of linear size L. The height of the interface at time t is represented by hi(t) and 
the number of bonds of each site by ni(t), with i — 1, 2, . . . , L . The initial condition 
is given by hi(0) = and ni(0) = 1, Vz, i.e., at the beginning of the simulation the 
interface is flat and a particle deposited into the site i will have ni = 1 bonds, Vi. A 
particle is deposited at random and search in its neighborhood for the site with the 
largest number of bonds and this site is chosen as the growing site. If the number of 
bonds of the particle cannot be increased in the neighborhood, the particle will choose 
the deposition site as the growing site and, if it has more than one site with the same 
number of bonds, one of them is chosen at random. One particle is incorporated to 
the growing site with a probability given by 



Pi = — (4) 
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where n max = 2d + 1 is the maximum number of bonds of one site, that is the 
case of a particle deposited in a site with its neighborhood fully occupied in a d- 
dimensional lattice. The variable v is the parameter of the model that control the 
growth probability of a given site and, for v = 0, the WV original growth model is 
recovered. Note that the rate growth of the sites with smaller neighbors number are 
decreased while the growth rate of sites with n — n max is always equal to 1 and an 
increase in the control parameter v emphasize the difference of these probabilities. 
After the deposition of a particle, the neighbors number of the growth site and of its 
neighbors must be updated. As usual |19) . we define one time unit as the tentative of 
deposition of L particles. 



3. Results 



As mentioned previously, if v = 0, the WV model is recovered as we can see in 
figures [lja) , for d = 1, and[5Ja), for d — 2. The interface obtained using this value 
exhibits a landscape governed by large plateaus separated by deep valleys that is the 
signature of models that incorporate diffusion [8]. For small values of v, for example 
v = 0.1 as shown in figure |TJb) , the interface resembles to that of [ija), however an 
attenuation of the largest steps are observed. Moreover, as can be seen from figure[TJc) 
and (d) and in figure [2jb) and (c) , an additional increase in the v value leads to an 
increase on this attenuation until the ridgeline landscape disappear completely, even 
for larger systems sizes and times asymptotically large. Therefore, in both (1 + 1) 




Figure 1. Interface morphology obtained after 10 s times step in a (1 + 1)- 
dimensions lattice of size L = 400. We used in these simulations v = 0.0, 0.1, 
0.50 and 2.0 for (a), (b), (c) and (d), respectively. 
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and (2 + l)-dimensions, the effect of the v parameter is to smooth the interface and 
this effect can be easily understood through an analyze of the growth probability. As 
observed previously, when we increase the v parameter, the rate growth of sites with 
n < n max decreases, hence the appearing of peaks are reduced. 

In order to characterize quantitatively the dependency of the roughening process 
by the parameter v, we consider the roughness of the interface that consist in the 
determination of the root mean square fluctuation in its height, i.e., 



!»(<) = 



i 



i N 

-J2Ht)-Ht)Y 



(5) 



where N — L d is the number of sites in the substrate, hi(t) is the height of the site i 





■L.,, 




(a) 



(b) 




(c) 



Figure 2. Interfaces obtained for distinct v values in a (2 + l)-dimensions system 
of linear size L = 100. We used v = 0.0 in (a), v = 2.0 in (b) and v = 4.0 in (c). 
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at time t and hit) is the mean height of the interface [Tj. We carried out simulations 
in systems of size ranging from L — 10 up to 10 5 for d = 1 while, for d = 2, we used 
systems from L = 10 up to 500. The v parameter was varied from to 8 and from 
to 4, for (1 + 1) and (2 + l)-dimensions, respectively. 

As we can observe in figure |3j for system in (1 + l)-dimensions, an important 
consequence of the introduction of the growth probability is the arising of different 
behaviors. If the v value is set to zero (upper curve in figure |3(a)T ), we find f3 w 0.375 
throughout all simulation, which means that the V 4 /i term of equation (J3J) dominates. 
However, for sufficiently small v values [v — 0.1), we observe a crossover in the /3 
exponent from a value ~ 0.375 to a value w 0.28, very close to j3 — 1/4 which 
means that the V 2 term of equation (|3| begins to prevail. This crossover is previously 
observed for the WV original model pj)l[2Q|. For a small increase in the v parameter, 
e.g. for v — 0.5 in figure 3(a)| the growth exponent is w 0.28 at all simulation. 
However, for an further increase in the v parameter (y — 2.0), we observe the exponent 
/3 w 0.33 after a crossover time. Finally, for largest v values (three bottom curve in 
figure |3~(b)[ ), we can observe (3 « 0.33 at all simulation. This new asymptotic growth 
exponent close to 1/3 is directly related to the rejection of that particle which choose 
a site with a very low growth probability. Therefore, this asymptotic value can be 
associated to the KPZ universality class, similar to the restricted solid-on-solid model 
proposed by Kim and Kosterlitz |19j . 

Corroborating the scaling of the roughness evolution, the saturation roughness 
also exhibits a crossover between different roughness exponents, as we can see in 
figure |4j When low values of the v parameter are used, for small values of L, the 
roughness exponent exhibits a value close to 3/2, as observed for the V 4 /i term, and 
a value close to 1/2 for largest values of the system size L. However, when large 
v values are used, this crossover disappear and the observed roughness exponent is 
always close to 1/2, a value obtained for both the EW (V 2 h term) and KPZ ((V/i) 2 
term) universality classes. 

In (2 + l)-dimensions we obtained a similar scene as that observed in (1 + 1)- 
dimensions. The figures [5] and [6] show the log-log and the semi log plots for the time 
evolution of the roughness and the system size dependency of the saturation roughness, 
respectively. As can be seen in these figures, for values of v < 0.5 a crossover from 
a power law to a logarithm growth was observed for both, the roughness and the 
saturation roughness scaling. The initial value obtained to the growth (/? 0.22) 
and roughness (a sa 0.78) exponents are very close to those found to the universality 
class of the V 4 h term, for wich (i — 1/5 and a = 2/3 and the asymptotic behavior is 
logarithmic (the growth and roughness exponent in the EW universality class are null 
that indicates a logarithmic behavior). For larger values of the scaling is always 
logarithmic. 

The results observed using intermediate v values for both, 1 + 1 and 2 + 1 
dimensions, indicating the EW asymptotic universality class are in good agreement 
with results found on literature [Till EH EH EH [21] ■ Besides, the time to the asymptotic 
behavior appear decreases as v is increased, so intermediate values of v can be used 
to anticipate the asymptotic universality class of the WV model. 



4. Conclusions 



In this work we studied a modified Wolf- Villain growth model with a probabilistic 
growth rule proportional to the power of the coordination number of a growing site, 
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i.e., p oc n v '. Simulations in (1 + l)-dimensions indicate that, for low v values, the 
model presents a crossover from the V 4 /i to the V 2 /i behavior, in agreement with the 
original model. We observe that the crossover time separating these two behavior 
decreases as v increases. This fact indicates that intermediary values of v can be 
used to anticipate the asymptotic universality class of the WV model. However, for 
intermediate v values, the model presents a crossover from the linear V 2 /i to the (V/i) 2 
nonlinear behavior and, for larger v values, the model falls into the KPZ universality 
class. The presence of the KPZ asymptotic universality class in this modified version 




Figure 3. Temporal roughness evolution for a (1 + l)-dimensions system of size 
L = 10 for various values of the parameter v as indicated in the legend. Notice 
that, to improve the visualization, the results for v < 2 arc show in (a), while for 
v > 2 are show in (b). 
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Figure 4. Saturation roughness as function of the system size obtained for a 
(1 + l)-dimensions system with distinct v values indicated in the legend of the 
figure. 




Figure 5. Roughness time evolution for distinct v values (0.0, 0.1, 0.5, 1.0 e 2.0) 
for a (2 + l)-dimensions system of linear size L = 500. In (a), a log-log plot is 
shown and, in (b), a semi-log plot. 



can be easily explained. The growth probability in this version introduce a refusal 
of particle. As know in the literature [TS], the refusal of particle is one of the basic 
condition to the appearance of the KPZ universality class. In (2 + l)-dimensions, the 
model presents the same type of crossover and, as in the (1 + l)-dimensional case, 
the crossover time is dependent of the v value. For larger value, we found that the 
model belongs to the EW universality class. In further work, we intend to apply the 
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Figure 6. Saturation roughness w sa t as a function of the system size L for 
distinct u values in a (2 + l)-dimcnsions system. In (a), a log-log plot is shown 
and, in (b), a semi-log plot. 



probabilistic rule to models, e.g. the DT model, that take into account the relaxation 
to the neighbor considering its coordination number. 
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